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G^ . Abstract 

T^ ' Results are presented for a combined analysis of the reactions pp — > pp, ird 






ird and vrd — > pp over the y/s interval from pion threshold to approximately 
2.4 GeV. These results for nd — > pp and ird elastic scattering are superior 

V 

\Q ' to our previous analyses of these reactions. In particular, the overall phase 

o' 

^^ [ in vrd — > pp has now been determined. Comparisons are made with previous 

1^^ I (separate and combined) analyses of this two-nucleon system. 
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I. INTRODUCTION 

An understanding of the A^A^ interaction is fundamental to studies of the more general 
nNN problem |]I[ . Below 1 GeV, in proton laboratory kinetic energy Tp for the A^A^ system, 
the dominant channels contributing to A^A^ inelasticity are nd and A^A [0]. At these energies, 
it is useful to employ a multi-channel formalism in analyzing all existing data simultaneously. 
In the present work, we have used the K-matrix formalism in order to unify the analysis 
of several reactions {pp ^ pp ^, nd -^ ird ^, and red ^ pp [0) which we have, in the 
past, considered separately. The range of y/s was chosen to include all of our results for the 
pion-induced reactions (T^r = — 500 MeV) . 

Clearly, we are not the first to consider this problem. A joint analysis of these three 
reactions, in a narrow energy range near the A^A threshold, was recently reported by Nagata 
et al. P|. This work used a mix of model-based and phenomenological results to investigate 
possible narrow structures in these reactions. An older work by Edwards [^ used the multi- 
channel K-matrix formalism to study the J^ = 2'^ and 3^ states associated with dibaryon 
candidates. 

The present analysis differs from those carried out previously in a number of important 
respects. We did not restrict our study to partial-waves containing interesting structures. 
For pp elastic scattering, all waves with J < 7 were used. Partial wave with J < 5 were 
retained for both red elastic scattering and Tcd -^ pp. In addition, the K-matrix parameters 
were determined solely from our fits to the available data bases for each separate reaction. 
No results of outside analyses or any model approaches were used as constraints. As a 
result, the amplitudes found in our K-matrix fits are as "unbiased" as those coming from 
the separate analyses [|]. 

In Section II, we will outline the K-matrix formalism used in this analysis. The combined 



and separate analyses will be compared in Section III. Conclusions and suggestions for 
further study will be given in Section IV. 

II. FORMALISM 



In order to analyze the reaction nd -^ pp along with elastic pp and nd scattering, we have 
constructed a K-matrix formalism having pp, nd and A^A channels. The energy-dependence 
of our global fit was obtained through a coupled-channel K-matrix form in order to ensure 
that unitarity would not be violated. The 'WA" channel is added to account for all channels 
other than pp and nd. The most important thresholds are illustrated schematically in Fig. 1. 
That this catch-all channel is indeed mainly A^A can be seen in Fig. 2, where the total cross 
sections for pp and nd scattering are broken into their components. 

As the elastic pp partial-wave analysis is far superior to the nd elastic and nd -^ pp 
analyses, we have carried out fits in which the pp partial-waves were held fixed. (The partial 
wave decomposition the of pp, nd, and A^A systems are given in Table I.) As described 
below, the pp amplitudes were used to fix some elements of the K-matrix, while the others 
were determined from a fit to the combined nd elastic and nd -^ pp data bases. 

States of a given total angular momentum and parity (J^) were parameterized by a 4x4 
K-matrix (Kj) which coupled to an appropriate A^A channel. Spin-mixed(2x2) pp states 
couple to unmixed nd states, and unmixed pp states couple to spin-mixed(2x2) nd states, 
so the nd — pp system is always represented by a 3x3 matrix. For example, the T-matrix 
{Tj) for J^ = 2+ (unmixed pp states) is given by 

pp nd_ nd^ 



( 'D2 'D2P 'D2F ^ 



To 
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pp 
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whereas the T-matrix for J^ = 2 (mixed pp states) is 

PP- pp^ nd 



To 



^P2 €2 ^P2D 


PP- 


£2 ^F2 ^F2D 


PP+ 
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The subscripts ± denote states with L = J ± 1. In the above, the mixing parameters 
(e) for elastic pp and ixd scattering are different. For the reaction Tid -^ pp, the notation 
{^Spp+^LJL^) of Ref. [0 is used. 

Adding an A^A channel results in a 4x4 T-matrix. Dropping the J-subscript, we write 
the K-matrix as 



/ 



K 
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\ 



(3) 



where Kpp is the elastic pp scattering sub-matrix, Kq and Kq are row and column vectors, 
and Ki is the sub-matrix of channels involving nd and A^A states. This K-matrix can be 
re-expressed as a T-matrix 



rp rp 

-Lpp Jo 



V 



To T, 



(4) 



using the relation T=K(l-iK) ^. We then have the correspondence 



T = K il —iK Y 
■'-pp ^^pp\^ ''^^ppj ) 



(5) 



where 



Kpp = Kpp + iKoil - iKi)-^ko 



(6) 



In order to ensure an exact fit to the pp elastic T-matrix, given by our most recent analysis 
of A^A^ elastic scattering to 1.6 GeV B, we take 



Kpp = Tpp{l + iTpp) - iKo(l - iKi) ^Kq. 



(7) 



The matrix elements are then expanded as polynomials in the pion energy times appropriate 
phase-space factors. The nd elastic and iid -^ pp T-matrix elements are extracted from Tq 
and Tj. 

III. PARTIAL- WAVE AMPLITUDES 

We have fitted the amplitudes for pp -^ pp and the existing data bases for nd -^ pp, 
and Tid — > vrd, using the K-matrix formalism outlined in Section II. The ird elastic and 
Tid -^ Tvd data bases used in this analysis are described in Refs. and 0, and available 
from the authors @]. The overall x^ fo^^ o^^ combined analysis is actually superior to that 
found in our single-reaction analyses. This is due to the improved parameterization scheme. 
A comparison is given in Table II. We should emphasize that the amplitudes for pp elastic 
scattering are the same as those given in Ref. p|. As mentioned above, this feature was 
built into our K-matrix parameterization. For this reason, we have omitted plots of the pp 
amplitudes |^. 

The results for nd elastic scattering are also qualitatively similar, up to the limit of our 
single-energy analyses. In Fig. 3 we compare the main partial-waves from our single-reaction 
analysis [^ and combined analysis (solution C500). Significant differences begin to appear 
above a pion laboratory kinetic energy of 300 MeV or 2.3 GeV in ^/s. (The ^D2 partial 
wave from C500 is an exception, departing from the single-reaction analysis near threshold.) 
The upper limit to our single-energy analyses is due to a sharp cutoff in the number of data. 
This is apparent in Fig. 2 of Ref. [^ . Much additional data above 300 MeV will be required 
before a stable solution to 500 MeV can be expected. 

A comparison of results for nd -^ pp reveals the most pronounced differences. One 
reason for this is the overall phase which was left undetermined in Ref. [0]. There, we 
arbitrarily chose the ^PiS wave to be real. In the present analysis, the overall phase has 
been determined. In Fig. 4 we show that the ^PiS phase is very different in the combined and 
separate analyses. Given the large difference in overall phase, we have chosen to compare the 



partial-wave amplitudes from the separate and combined analyses in terms of their moduli. 
This comparison is made in Fig. 5. As was the case for nd elastic scattering, differences are 
most significant above approximately 2.3 GeV in ^/s. A similar lack of data exists above 
this energy. 

In general we see a good agreement for the dominant amplitudes found in the separate 
and combined analyses. Figures 3 and 5 also display our single-energy analyses which were 
done in order to search for structure which may be missing from the energy- dependent 
fit. (Details of the single-energy analyses are given in Refs. [|],0|.) A comparison of the 
single-energy and energy- dependent fits is given in Tables III and IV. 

IV. SUMMARY AND CONCLUSIONS 

We have obtained new partial- wave amplitudes for ird elastic scattering and the reaction 
7id —>■ pp, using a K-matrix method which utilized information from our elastic pp scattering 
analysis. In addition to producing amplitudes more tightly constrained by unitarity, we have 
resolved the overall phase ambiguity existing in our previous ird —>■ pp analysis. 

As mentioned in Section III, the combined analysis has resulted in a slightly improved 
fit to the nd elastic and nd -^ pp data bases. The most noticeable differences, at the partial- 
wave level, appear at higher energies where the existing data are sparse. It is difficult to 
find cases where the fit has been dramatically improved. One exception is the set of ird 
total cross section data between 300 and 500 MeV. Here the combined analysis is much 
more successful in reproducing the energy dependence. The combined analysis gives total 
cross sections which begin to rise at 500 MeV, whereas the separate analysis shows a fairly 
monotonic decrease from 400 to 500 MeV. The behavior seen in the combined analysis seems 
reasonable, as the nd total cross sections do begin to rise just beyond the upper energy limit 
of our analysis. Many of the individual partial-wave amplitudes from C500 show rising 
imaginary parts near 500 MeV, a feature absent in the analysis of nd elastic data alone. 

The present analysis has also resulted in a unified description of the resonancelike behav- 



ior previously noted in our separate analyses of pp and 7rd ^ elastic scattering, and the 
reaction nd —>■ pp ^. This behavior [|I^ has been variously described as "resonant" (due to 
the creation of dibaryon resonances) and "pseudo-resonant" (due to the A^A intermediate 
state). We expect that our combined analysis will further constrain models based on these 
two mechanisms. 
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Figure captions 

Figure 1. Energy scale in terms of the total center-of-mass energy (i/i) and the inci- 
dent kinetic energies of the pp (Tp) and nd (T^r) initial states. The locations of 
relevant thresholds are also displayed. 

Figure 2. Total cross sections atot (solid) and total elastic cross sections a^i (dashed) 
correspond to the C500 solution. Data for atot (open circles) are taken from the 
SAID database [Q. (a) Dash-dotted lines, corresponding to the C500 solution, 
show the total cross sections (aTrd) for pp —>■ nd. The corresponding data from 
the SAID database [Q are plotted as open triangles. The remainder (Aa) is 
given by atot — osi — CTnd and plotted as a dotted hne. Total cross sections for 
the reactions pp -^ A^p + A++n are plotted as dark circles, (b) Dash-dotted 
lines (C500) show the total cross sections {app) for ird —>■ pp. The corresponding 
data from the SAID database 0] are plotted as open triangles. The remainder 
(Act) is given by atot — (^ei — (^pp and plotted as a dotted line. 

Figure 3. Partial-wave amplitudes of the reaction nd —^ nd from T^ = to 500 MeV. 
Solid (dashed) curves give the real (imaginary) parts of amplitudes corresponding 
to the C500 solution. Our previous analysis (SM94) p is plotted with long dash- 
dotted (real part) and short dash-dotted (imaginary part) lines. The dotted 
curve gives the value of Im T - T^ - T^r, where T^/- is the spin- flip amplitude for 
C500. The real (imaginary) parts of single-energy solutions are plotted as filled 
(open) circles. All amplitudes have been multiplied by a factor of 10'^ and are 
dimensionless. Plotted are the dominant partial-wave amplitudes: (a) '^Pq (0^), 
(b) ^S, (1-), (c) ^P, (2+), (d) ^D, (2-), (e) -'D, (3"). 

Figure 4. Comparison of the ^PiS partial waves for nd —>■ pp obtained in the separate 
and combined fits. The real (imaginary) part of solution C500 is plotted as 

10 



a solid (dashed) line. The purely real partial wave from our separate analysis 
(SP96) [0 is plotted as a dot-dashed line. 

Figure 5. Moduli of the partial-wave amplitudes for nd -^ pp from Ttt = to 500 MeV. 
The solid and dashed curves give the amplitudes corresponding to the C500 and 
SP96 10] solutions respectively. Moduli of the single-energy solutions are plotted 
as filled circles. All amplitudes have been multiplied by a factor of 10^ and are 
dimensionless. Only dominant partial- waves have been plotted: (a) ^SqP (0^), 
(b) ^P,S (1~), (c) 'D^P (2+), (d) ^P^D (2-), (e) 'F^D {3'). 
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TABLES 
TABLE L Partial wave decomposition of pp, ird, and A'^A systems. 



JP 


vrd 


pp 


NA 


0+ 


'Po 


'So 


'Do 


0- 




'Po 


'Po 



3Pi 3^1, ^Di 

^Si, ^Di ^Pi ^Pi, ^Fi 

3P2, ^F2 ^D2 3^2 

^P2, ^F2 'D2 5^2, 'D2 

^P2, ^F2 'D2 'D2, 5^2 

3n 3p3p 3p3p 

^2 ^2, ^2 ^2, ^2 

''D2 ^P2, ^F2 ^P2, 'F2 

'F, 3^3, 'Gs 

'F3 ^^3, 'G3 

'i^3, 'Gs 'i^3 'i^3, ^Fs 

3 7-) 3/^ 3/7 5 p 5 771 

^3, Lr3 -f^3 -rs, ^3 

3 7-) 3/^ 3 77 5 77 5 rr 

-t>'35 L^3 ^3 -f^3, lis 

3^4, 37^4 ^G4 3^4 

3^4, 3/^4 IG4 ^A, ^G4 

r^, X34 Gr4 ^4] -'4 

3/^ 3773f7 377317 

Lt4 -r4, -714 -^4) -"4 

3G4 ^^4, 3^4 ^F4, ^i^4 
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TABLE II. Comparison of the combined analysis (C500) and our previous (separate) analyses. 
WI96 for pp ^ pp [||, SM94 for vrd ^ vrd |§, and SP96 for vrd -^ pp 0. The relevant energy 
ranges are: T^ = 0-500 MeV, Tp = 288-1290 MeV, and ^/s = 2015-2440 MeV. 



Reaction 



Separate 
X^/Data 



Combined 
X^/Data 



pp -^ pp 
vrd — > vrd 
nd — > pp 



17380/10496 
2745/1362 
7716/4787 



17380/10496 
2418/1362 

7570/4787 
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TABLE III. Comparison of single-energy (binned) and energy-dependent combined analyses of 
TTcl elastic scattering data. Nprm is the number of parameters varied in the single-energy fits. Xe 
is due to the energy-dependent fit (C500) taken over the same energy interval. 

T^ (MeV) Range (MeV) Nprm xV^ata x| 



65 


58.0 - 72.0 


2 


106/54 


102 


87 


72.0 - 92.0 


6 


20/24 


21 


111 


107.5 - 125.2 


10 


68/82 


66 


125 


115.0-134.0 


12 


155/170 


184 


134 


124.0 - 142.8 


14 


315/258 


344 


142 


133.0-152.0 


16 


356/284 


397 


151 


141.0-160.6 


16 


193/154 


216 


182 


174.0 - 189.5 


18 


302/168 


396 


216 


206.0 - 220.0 


18 


158/99 


200 


230 


220.0 - 238.0 


18 


64/53 


111 


256 


254.0 - 260.0 


16 


132/125 


185 


275 


270.5 - 284.4 


16 


22/40 


42 


294 


284.4 - 300.0 


16 


267/132 


324 
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TABLE IV. Comparison of single-energy (binned) and energy-dependent combined analyses of 
TTcl -^ pp reaction data. Nprm is the number of parameters varied in the single-energy fits, xl; is 
due to the energy-dependent fit (C500) taken over the same energy interval. 

T^ (MeV) Range (MeV) Nprm xV^ata x\ 



25 


12.8-37.4 


10 


527/241 


542 


50 


37.6 - 60.7 


12 


188/168 


205 


75 


62.9 - 87.3 


14 


590/426 


628 


100 


91.0- 114.0 


14 


1263/611 


1379 


125 


113.8-137.1 


16 


729/512 


756 


150 


140.0 - 162.0 


20 


743/630 


792 


175 


165.0 - 187.3 


22 


343/280 


426 


200 


191.3-210.3 


20 


120/193 


153 


225 


217.9 - 235.9 


22 


217/229 


291 


250 


238.9 - 262.0 


22 


595/483 


685 


275 


264.9 - 285.1 


22 


204/109 


280 


300 


291.6 - 307.4 


24 


198/212 


235 


325 


318.9-330.0 


24 


142/161 


234 


350 


341.4 - 360.3 


24 


201/185 


233 


375 


371.4 - 375.7 


24 


32/26 


42 


400 


390.0 - 400.0 


24 


19/28 


34 


425 


417.0 - 420.0 


24 


50/28 


55 


450 


437.6 - 456.5 


22 


122/48 


231 


475 


473.8 - 487.4 


22 


24/24 


39 


500 


495.9 - 506.5 


22 


49/45 


281 
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